Two exact properties of the perturbative expansion for the two-dimensional electron 
liquid with Rashba or Dresselhaus spin-orbit coupling 
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We have identified two useful exact properties of the perturbative expansion for the case of a two- 
dimensional electron liquid with Rashba or Dresselhaus spin-orbit interaction and in the absence of 
magnetic field. The results allow us to draw interesting conclusions regarding the dependence of the 
exchange and correlation energy and of the quasiparticle properties on the strength of the spin-orbit 
coupling which are valid to all orders in the electron-electron interaction. 
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I. INTRODUCTION 

Modern yearning for understanding the electronic 
properties of devices made out of materials in which 
spins play an important role has rekindled interest in 
the study of the two-dimensional electron liquid in the 
presence of Rashba or Dresselhaus spin-orbit coupling*^— 
To date, the effects of the electron-electron interac- 
tions in these system have been studied by means 
of approximated methods, notably the random-phase- 
approximation (RPA) and the mean-field theory^ - — 
Quasiparticles parameters were studied in Rcfs. while 
Refs. @4l3 focused on the ground-state properties and 
mean-field phase diagram. Recent studies of spatially in- 
homogeneous Overhauser's type of instabilities (chiral- 
ity, charge and spin) and of the spin-susceptib ility in 
these systems can be found in Rcfs. [ll| and Ibi - fla re- 
spectively (including non-analytic corrections 15 ). Cal- 
culations of realistic systems of lower dimensionality in 
the presence of spin-orbit coupling have also been car- 
ried out by means of numerical density functional theory 
methods! 1 : 18 However, for want of more accurate ex- 
change and correlation functionals, the latter have been 
implemented by employing formulas for the energy of the 
electron liquid obtained in the absence of spin-orbit cou- 
pling, an approximation that a priori could appear at 
best rather crude. Only recently, quantum Monte-Carlo 
numerical data of the energy were obtained for electrons 
in two dimensions with Rashba spin-orbit coupling^ 

In this paper we analyze some formal aspects of the 
perturbative expansion in the electron-electron interac- 
tion, without magnetic field and under the assumption 
that the system will behave as a Fermi liquid. In doing 
so we will show that certain exact identities lead to the 
conclusion that, at least in the high density regime or 
for small spin-orbit coupling strength, the effects of the 
interactions are only marginally influenced by the spin- 
orbit coupling. In particular we will show that, to all 
order in the electron-electron interaction, the corrections 
to the total energy and the Fermi liquid properties for- 
mally vanish to lowest order in the spin-orbit coupling. 



As consequence, our results lend some degree of justifi- 
cation and theoretical underpinning for many of the nu- 
merical studies carried out by neglecting such effects. 

In Ref. [IH a canonical transformation was proposed for 
the case of small dots that allows one to renormalize the 
hamiltonian into one in which a generic combination of 
Rashba and Dresselhaus spin-orbit coupling terms only 
appears in second order. The vanishing of such spin-orbit 
interaction effects at small coupling was later verified 
within the RPA in Ref. |5| for the self-energy. Our analysis 
not only extends these results to the fully interacting case 
but also explicitly shows how the repopulation of momen- 
tum space with respect to the non interacting situation- 
is responsible for the renormalization of the Rashba and 
Dresselhaus coupling. Our formalism also sets the stage 
for a systematic analysis of the perturbation theory to all 
order for regimes in which the series is expected to con- 
verge. The specific application to the case of Coulomb 
interaction in the high-density regime and a comparison 
of analytic formulas for the exchange and correlation en- 
ergy to the numerical Monte-Carlo results of Ref. [ID can 
be found in Ref. [T^- 

The paper is organized as follows: In Section [TTJ we 
provide the definition of the model and the necessary no- 
tation; In Section Hill we prove an exact identity concern- 
ing the dependence of the total interacting energy on the 
suitably defined density dependent bare spin-orbit cou- 
pling and the chirality, the relevant Fermi surface repop- 
ulation parameter; In Section HVl a similar relationship is 
obtained for the self-energy perturbative series; Finally, 
Section IVl provides a discussion of our results. 



II. GENERALITIES AND BASIC FORMALISM 

Consider the many-body problem associated with the 
hamiltonian 
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where the presence of an homogeneous neutralizing back- 
ground is understood. The second term in the square 
brackets represents a spin-orbit coupling of the linear 
Rashba type and can be readily seen to be equivalent 
to a Dresselhaus spin-orbit coupling term»i~— While in 
Eq. (p} the interaction is Coulombic, as it will be shown, 
our results are valid for a generic interparticle potential. 
The single particle solutions of the non interacting sys- 
tem are the well known chiral states: 
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FIG. 1. Left: interaction vertex, associated with a (k/x|k'//) 
factor. Right: example of a fermionic loop. 



where L is the linear size of the system, </>k is the angle 
between the wave vector k and the x axis and we have 
defined the spinors |k±). 

It is useful at this point to simplify the notation by 
rescaling all wave vectors by /cp = v2~7rn, with n the areal 



density, and all energies (frequencies) by ep 
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the 



Fermi energy in the absence of spin-orbit coupling ( ^& ) . 
With the present notation the eigenenergies correspond- 
ing to © can be simply expressed as 



e±(k) = k 2 Tgk , 
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where we have defined the dimensionless coupling con- 
stant 
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which plays an important role in the following. Without 
loss of generality, we assume g > 0. 

As discussed in Ref. the unperturbed, spatially uni- 
form states with symmetric k-space occupation can be 
completely characterized by the generalized chirality x- 
This parameter is defined by the Fermi surfaces, two cir- 
cles with radii given by 



fc± = y/\T± 



x\ 



(5) 



When both chiral bands are occupied this quantity coin- 
cides with the regular chirality: 
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where N± is the number of electrons in each chiral band. 
This is the case for not too low densities and it is the 
situation we will deal with here.- 

The corresponding occupation numbers can then be 
expressed as follows: 



n±(fc) = 0(yi±x-fc) , forfj<x<l, (7) 

where 6(x) is the usual step function. The non interact- 
ing energy per particle (in units of ep) is in turn given 
I'v 



Here the first term is the kinetic energy while the sec- 
ond one is the spin-orbit coupling contribution^ For 
non interacting electrons the value of x that minimizes 
this expression depends only on the strength of the spin- 
orbit coupling and the electron density and, remarkably, 
is in fact uniquely determined by the parameter g [de- 
fined above in Eq. Q] as follows: 
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for < g < y/2 
for g > V2 . 
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The perturbative analysis is very similar to the stan- 
dard case without spin-orbit coupling.— The fully in- 
teracting Green's function G fl (k,t) = — i{Tl^^(t)b^ {0)) 
(fi = ±) is obtained as an average of the time ordered 
product of creation (destruction) operators SL, (&k^) of 
the chiral states of Eq. ^ The only modification to the 
diagrammatic rules is that each vertex, beside the delta 
functions ensuring frequency and momentum conserva- 
tion, is associated with the scalar product (k/z|k'//) [the 
spinors having been defined in Eq. fl3J)]. This is exem- 
plified in Fig. [TJ Accordingly, the overall dependence of 
the diagram on the variables x an d g is solely contained 
in the non interacting Green's functions. The Fourier 
transform of the latter is given by 



G 0fl (k,uj) = 
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where n±(k) is defined in ([7]) and 77 = + is an 
infinitesimal quantity. Notice that the dependence 
on x is fully contained in the occupation numbers, 
while the dependence on g is only contained in e M (fc). 
A similar set of diagrammatic rules applies to the 
self-energy, defined through the relation G M (fc,cj) = 
l/(w-e M (fc) + £ M (fc,u;)//i). 



III. AN EXACT PROPERTY OF THE ENERGY 
EXPANSION 
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For any given value of the chirality, the fully interacting 
total energy per particle can be obtained from the follow- 
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ing integration over the coupling constant formula^ 
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where the non interacting contribution was given in 
Eq. ((H) . In the second term the task is to integrate over A 
a sum of closed diagrams, each multiplied by a factor of 
A" -1 , where n is the order of the interaction. Consider 
now the expression D of one of these closed diagrams, 
contributing to the second term of Eq. (fTTj) . It is the 
purpose of this Section to show that 



(12) 



d 2 D 




d 2 D 


8 2 D 


d 2 x 





d 2 g 


o d xdg 



where the derivatives are evaluated at x — .9 = 0. Our 
treatment is carried out to all orders of perturbation the- 
ory, for a generic diagram D and a general two-body po- 
tential v(q). 

To prove Eq. (fT2|) . we begin by noticing that all 
the closed diagrams are in general comprised of several 
fermionic loops. Let us focus our attention on the con- 
tribution to Eq. (fl2|) stemming from of one such loop 
(schematically illustrated in Fig. [TJ. If the loop contains 
N solid lines we can write: 



D 



/ME/^/£rK<^.> 

{Mi} i=1 



x (pifiilpi+im+i) , (13) 

where {/ii} = {fix, f^2, ■ • • Miv} and N + 1 = 1. As in 
Fig. [TJ the internal momenta and frequencies are pi = 
k + qi and = co + ili where k and ui are the momentum 
and the frequency flowing in the loop. Finally the square 
bracket [. . .] represents the remaining expression of the 
diagram, which is independent of k and u and does not 
enter the derivation below. 

We proceed by examining the first derivatives of 
Eq. (IT3l) with respect to the variables x an d g for g — 
X = 0. It is readily seen that D depends on x and g only 
through the non interacting Green's functions (flU)) . As 
a consequence, a direct inspection of Eq. (13) leads to 
the conclusion that both first derivatives for g = x — 
are given by a sum of N terms, of which the j-th term 
contains the first derivative of the corresponding Green's 
function in the diagram. Since both 
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are proportional to /ij, this yields the factor 
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which, as indicated, vanishes upon spin summation. 

The second derivatives of D have more complicated 
expressions that are in general non-vanishing also for 
X = g = 0. We observe next that the dependence of 



Go^ on x an d g is such that -g- only acts on the occu- 
pation numbers n±(k), while ^ only acts on the energy 
denominators of this function. To prove Eq. 
then introduce as a formal device the quantity 



we 



D„ 



duj dk x dk g 
2^ (2tt) 2 



S(k x 



N 



(15) 



where each of the the indexes 



and 



= k„ 



i=X 

b can be % or j while 
p Xi i = k x + qj and p g ^ — k g + q ; . This expression 
corresponds to the same diagrammatic contribution of 
Eq. (|13p . The square parenthesis in the first line co- 
incides with the expression omitted in Eq. (|13[) . but is 
here evaluated at g = x = 0. We have also defined the 
'resolved' non-interacting Green's function for vanishing 
spin-orbit coupling (g = x — 0), i.e. 



g (k x ,k g ,u) = 



l-n (k x ) n (k x ) 



to — k 2 + ir\ to — k 2 — irj 
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where no(k) = 9(1 — k). In this expression, the momenta 
fc x and kg respectively appearing in the occupation num- 
bers and the energy denominators are treated as two in- 
dependent variables. Notice that, except for the differ- 
ential operator inside the square parenthesis, Eq. f| 15[) 
coincides with the X = .9 = expression of the diagram 
D. At this point, by integrating by parts (so that the 
derivatives act on the two dimensional delta function), 
one can see that 



Dxx ~ Dgg — D xg 



(17) 



The desired identity Eq. (Ti"2"j) follows from the relation 
d 2 D 



dadb 



Dab 
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which is proved next. 

To derive Eq. (fT8|) we consider first the equivalence of 



d 2 D 
dx' 2 



with D y 



' Xx . The former is the sum of several terms, 

of which some contain the second derivative of a single 
Green's function, and others contain the product of first 
derivatives of two distinct Green's functions. The same 
holds true for D xx [see Eq. (|15j) ]. The euqlity between 

corresponding second derivative terms of ^§ and 



can be established by using the relation 

_ d 2 G 0il {p,u) 



o 



G (p X ,Pg,u) 
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dx 2 



(19) 



where on the left side we set p x = p g = p and on the right 
side g — X — 0- The necessary factor of 2 is obtained 
performing the spin summation 
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For the first derivative terms, we can use the formula 

(« = 9,X) 



g (p x ,p g ,uj) 



dG 0fl (p,uj) 



da 



Px=Pg=P 
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P 
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The factors P J f k are recovered by performing the follow- 
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ing spin summations 

N 



This establishes the desired result 0^ 



P 3 ■ Pk 

P 3 Pk 

= D 



(22) 



The equality 



d 2 D 



dxdg 



D xg can be proved in a similar 



way. For this case, the following relation proves useful 
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Finally, to prove the equivalence of 
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use the following formula 
GoiPxiPg^) 



d 2 G ^(p,uj) 



Px=Pg=P 



dg 2 



and D gg , we 



-[00 (p, w)] 



(24) 

where Go(p, u) is the usual non interacting Green's func- 
tion without spin-orbit coupling (g = \ = 0)«22, A slight 
complication arises, due to the [t/o(Pjw)] 2 term of (|24| . 
This produces an additional contribution of the following 
form 



dk 

(2^0 



However, (|25p is seen to vanish upon integration over the 
loop frequency tu. 



IV. AN EXACT PROPERTY OF THE 
SELF-ENERGY EXPANSION 

The fully interacting self-energy Yt^kjUj) in the pres- 
ence of spin-orbit coupling satisfies a similar exact re- 
lationship to linear order in g. Again, we consider its 
diagrammatic expansion and show that for any given di- 
agram £> s : 



fi dDT, (k,uj) 
2 dk 



(26) 



where Eoffc, ui) is the interacting self energy without spin- 
orbit coupling (therefore, D^ is obtained setting g = 



X = in Ds^). We will derive Eq. ([26]) assuming that 
X — g and, as in Eq. (|T2l . the left side derivative is 
evaluated at g = 0. Within the notation of the previous 
Section, D^, is written as 



N 

En^o 



JY 



'Owtp^i) Yl(PjH\Pj+iH+i) ■ 

3=0 

(27) 

Here, only Green's functions connected which the exter- 
nal momentum k are explicitly written. If 1 < i < N we 
set pi = k + q, while p = pjv+i = k and fi = fijsr+i = 
(i. 

Also in this case the derivative -§- involves all the non 
interacting Green's functions appearing in the diagram. 
On the other hand, as we have discussed in the previous 
Section, one needs not to worry about the (vanishing) 
contribution of all the sets of Green's functions that are 
involved in fermionic loops. 



Therefore, the expression of 



dg 



consists of the sum 



of N terms, involving the first derivatives of the Green's 
functions explicitly appearing in Eq. (|27[) . The final re- 
sult is 
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where the factor Qo(j)j,ojj) is omitted in the second line. 
We also notice that, in the first line, the expression in 

the round brackets is independent of fij (since 



i)G„ 



fi). Furthermore the sums over the spin indexes {fJ-i} 
{[11,11,2, . . . [In} were performed using the relation 
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We then consider the expression Ds of the same di- 
agram occurring in the expansion of the interacting self 
energy in the absence of spin-orbit coupling: 



N 



;.=i 



Its derivative can be directly compared to Eq. 
making use of the relation 



dk 



dGo^ipj,^) 



which is easily verified for \ 
ately follows. 
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V. DISCUSSION 

We have derived by diagrammatic means two exact 
relationships involving the total energy expansion and 
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the self-energy of a two dimensional system in the pres- 
ence of Rashba or Dresselhaus spin-orbit coupling. The 
first identity is given in Eq. (fT2|) and involves the second 
derivatives of the energy diagrams with respect to the 
density dependent dimensionless spin-orbit coupling g, 
defined in Eq. (j4} , and the chirality x, defined in Eq. ((6|) . 
Notice that, as made clear in Eq. ([S]), \ takes into ac- 
count the renormalized momentum space repopulation 
brought about by interactions 4 Since Eq. (fl"2j) is valid 
for any diagrammatic contribution, we can formally sum 
the perturbative expansion in the electron interaction. 
This leads to the following formula, for the small g and 
X expression of the extra contribution to the electronic 
exchange and the correlation energy due to the spin-orbit 
coupling: 

5£ xc (9,x)=C(g- X f + -- - , (32) 

where the constant C is unspecified, and terms of higher 
order in g and x are omitted. It is clear from this result 
that the total energy is an extremum for x — 9- Com- 
paring with Eq. (0) we conclude that to linear order in 
g the interactions lead to no change in i.e. no repop- 
ulation, and there is no energy correction S£ xc (g,x) to 
quadratic order in g. 10 This remarkable result does not 
hold for other types of spin-orbit coupling^ 

Our second result, Eq. (|26| . involves the self-energy 
and allows us to write that to linear order in g 

E M (A,a;) = So(A!-^ J w) + ... . (33) 

From this relation we conclude that, to linear order in 
g, all the quasiparticles properties (effective mass, life- 
time, . . . ) on the Fermi surfaces k± ~ 1 ± | can be 
simply obtained from the case without spin-orbit cou- 
pling. Equation (|3"3"]l extends the results of Ref. d where 
this statement was derived in the RPA approximation 
for Coulomb interactions. We observe that, similarly to 
the case of the phonon self-energy in a Fermi liquid, the 
interaction self-energy for the case of spin-orbit coupling 
is seen to "ride" the Fermi surface. 



Finally, we would like to comment about the relation 
of our calculation to non-analyticities, which are known 
for a long time to occur in the theory of the electron 
liquid 24 and have attracted much attention in more 
recent years (see, for example, Refs. Ila. 12614311 ). Non- 
analytic terms in the interactions are fully included in 
our results, since our derivation applies to all orders in 
perturbation theory. However, an analytic expansion in 
the spin-orbit coupling parameters, g and x> is assumed. 

This procedure is justified here because we restrict our- 
selves to low orders in the spin-orbit coupling parameters 
(first-order for the self-energy and second order for the 
total energy). For the total energy, non-analytic correc- 
tions can only appear starting from third order in the rel- 
evant energy scale, as it can be concluded from a power- 
counting argument i 30 ' 31 For example, a contribution to 
the total energy oc \o\E\ (where Ez is the Zeeman en- 
ergy) is responsible of a non-analytic correction cx |a| 
to the spin-susceptibility^ This contribution to the en- 
ergy is vanishing in the absence of magnetic field but 
a term oc |a| 3 ln 2 \a\ can be derived at fourth order in 
the electron interaction^ A non-analytic correction to 
the exchange energy is also obtained in Ref. UH, propor- 
tional to x 4 ln|x|. For the self-energy, no terms linear 
in a were found in Refs. [H and [13 and non-analytic cor- 
rections proportional to a 2 In |a| were derived in Ref. [!. 
All these non-analyticities do not affect our main results, 
Eqs. (|3^|) and (|33p . since they only appear at higher order 
in the spin-orbit coupling. 
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